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Abstract 

In samples used to maintain the US resistance standard the breakdown of the dissipationless integer quantum Hall effect 
occurs as a series of dissipative voltage steps. A mechanism for this type of breakdown is proposed, based on the generation 
of magneto-excitons when the quantum Hall fluid flows past an ionised impurity above a critical velocity. The calculated 
generation rate gives a voltage step height in good agreement with measurements. We also compare this model to a 
hydrodynamic description of breakdown. 
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In the integer quantum Hall effect (QHE) regime, a 
two-dimensional electron fluid carries an almost dissi- 
pationless current and the ratio of the current, to the 
Hall voltage is quantized in units of e 2 /h. However, 
above a critical current, the dissipative voltage, Vx, 
measured along the direction of current flow, increases 
rapidly, leading to breakdown of the QHE. For certain 
samples, including those used to maintain the US resis- 
tance standard at the National Institute of Standards 
and Technology (NIST) [1,2], breakdown occurs as a 
series of steps, in V x , of height AV X ~ 5mV . Below we 
examine this step height using two distinct methods. 
First we perform a fully quantum calculation for the 
rate of formation of e-h pairs due to a single charged 
impurity, as a function of electric field, E y . In doing 
this we obtain a relationship between V x and E y . We 
find that as E y increases V x increases rapidly to a maxi- 
mum and then falls slowly. We then argue that the only 
stable points on this graph are firstly when V x is less 



than the small dissipative voltage governed by o~ xx and 
secondly when V x is a maximum. Hence we are able to 
estimate the step height AV X . We then consider a hy- 
drodynamical model [3] of this process starting from an 
effective Euler equation. By examining small pertuba- 
tions in the velocity field we find a direct analogy to the 
elastic inter-Landau level tunneling condition [4] , from 
which we can again estimate AV X . Finally, we compare 
the results of the two approaches with experiment and 
comment on how this work relates to bootstrap elec- 
tron heating (BSEH) [5] model of the breakdown of the 
integer QHE. 

The starting point is the Fermi golden rule, which we 
use to calculate the rate of generation of e-h pairs, W , 
due to a single charged impurity in a uniform electric 
field. Using the Landau gauge, this gives 
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Fig. 1. The magneto-exciton mode energy, 
fojj c + e 2 /(4yTKis)Ao i (i) [Q] (solid line) for excitations from 
ra = to n = 1, for E y = 0; the dashed line is the electro- 
static energy, eE y l 2 B Q. B = 12. 3T and E y = 1.5 X 10 6 V/m. 
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where V(r) is the impurity Coulomb potential, 
4>n(r±,k x ) is the electronic eigenfunction in the xy 
plane and Ho (2) is the envelope function of the first 
electronic sub-band. We assume that the n lower Lan- 
dau levels are filled and the (n + 1) level is empty. 
Hence there is no static screening of the impurity 
charge. Calculating the transition rate out of state 
(n,k x ), by summing over k' x in Eq. (1) we obtain the 
following energy conservation condition 
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Ae = hu>c+-. — A n (n+1) [k x -k x ]-eE y l%(k' x ~k x ) = 0, 
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where uj c is the cyclotron frequency, h is the magnetic 
length, « is the dielectric constant, E y is the y com- 
ponent of the Hall electric field, h(k' x — k x ) = HQ is 
the momentum change and A Tlj ( n+1 ) [Q] includes the 
local and exchange field corrections [6] for an excita- 
tion from (n, k x ) to (n + l,k' x ). The condition given 
by Eq. (2) is more precise than that presented in Ref. 
[4] for elastic inter-Landau level tunneling. The im- 
provement comes from inclusion of local and exchange 
field corrections for the excitation process, i.e. in Ref. 
[4] A n ( n+1 )[Q] = 0. When the dashed line crosses the 
solid curve, in Fig. 1, the condition given by Eq. (2) is 
met. At this point, energy is conserved in moving an 
electron from the filled lower Landau level (LL) to the 
unoccupied upper LL. For a given electric field, E y , we 
then use Eq. (1) to calculate the rate of generation of 
e-h pairs due to the presence of an impurity. 

Consider a local region of the sample where E y is 
large enough to create e-h pairs at a given rate, due 




Fig. 2. A schematic diagram showing the variation of the 
Landau level energy across the Hall bar. The e-h pairs are 
formed in region I by scattering from a charged impurity 
at a rate W . The cross-section of the bar is split into three 
regions. In regions I and III, near the edges of the bar, the 
voltage drop is high compared to region II in the centre of 
the hall bar. 

to scattering from a charged impurity. A pair created 
close to an impurity drifts along the Hall bar at a veloc- 
ity E y / B, so for a fixed generation rate, a stream of e-h 
pairs moves along the Hall bar. Then, for a filling fac- 
tor of 2 we have a situation in which a small fraction of 
electrons in the lower LL (n — 0) have been replaced by 
holes and the previously empty upper LL (n = 1) con- 
tains some electrons. As the e-h pairs move away from 
the high field region, the spacing between the electron 
and hole in a pair increases and most pairs eventually 
ionise by acoustic phonon emission. Due to the absence 
of empty states into which the excited electron can re- 
lax and neglecting weak, second order Auger processes, 
we can assume that all the generated e-h pairs even- 
tually ionise, leading to a dissipative current i = eW, 
where W is calculated from Eq. (1). This process is 
shown schematically in Fig. 2. From this figure we see 
that an e-h pair created close to the left (hot) edge of 
the Hall bar leaves a hole which ionises at the hot edge 
of the bar whilst the excited electron, in the n + 1 Lan- 
dau level, relaxes its energy due to phonon emission as 
it passes through regions I, II and III. At filling factor 
v — 2, this process gives rise to a dissipative voltage 



The NIST experiments carried out at v = 2 and 
B — 12.3T show a series of up to 20 dissipative steps 
in V x of regular height AV X ~ 5mV. In Fig. 3 we 
plot our calculated dissipative voltage versus the back- 



2 



ground electric field E y . Since W is strongly influenced 
by the overlap between the wave-functions in the occu- 
pied (n = 0) and unoccupied (n = 1) LLs, it increases 
rapidly at a critical electric field. This occurs when V x 
is comparable to the small background dissipative volt- 
age governed by a xx . The rate, W, and hence V x , then 
reaches a maximum when the electric field is such that 
the e-h pairs are formed close to the roton minimum of 
the magneto-exciton dispersion curve. For a given E y , 
for which W is finite, the pairs formed at the break- 
down point relax, so that electrons in the upper LL 
moving towards one side of the Hall bar, whilst holes 
move in the opposite direction. The presence of these 
carriers screens the Hall field over much of the Hall 
bar. Since the Hall voltage in the NIST experiments 
remains constant at its quantized value over the mag- 
netic field range where the dissipative steps occur, this 
screening effect tends to enhance the electric field at 
the breakdown point. Thus, as the critical electric field 
is reached, the generation rate at the breakdown point 
increases rapidly, inducing a further increase in E y due 
to the screening of the Hall field in other regions of 
the Hall bar. We can see how this happens in Fig. 2. 
The electron is excited in region I and then decays due 
to phonon emission to region II. The voltage drop, V2, 
is small so the electron takes a relatively long time to 
traverse this region. Hence, if the excited electrons in 
region I are generated at a constant rate, W , they tend 
to accumulate in region II. Since Vh remains constant 
in the magnetic field range over which the steps are ob- 
served [2] and the accumulated electrons in region II 
screen the local electric field in this region the voltage 
drop, V2, must reduce and hence the voltage dropped 
in regions I and III must increase. As a result, the elec- 
tric field in regions I and III increases. This increase 
in E y increases W until it reaches a maximum. Thus 
for breakdown at a single charged impurity the system 
should switch between two stable states, correspond- 
ing to V x = and V x = 5.6mV, which corresponds 
to the maximum value of V x in Fig. 3. This is in good 
agreement with NIST experiments. In the experiments 
a series of up to twenty steps are observed and we at- 
tribute each step to the formation of separate streams 
of e-h pairs generated by different impurities. 

We now examine more closely an analogy recently 
proposed by one of us [3] between the breakdown of the 
QHE and the breakdown of laminar fluid flow around 
an obstacle. Our starting point for such a model is an 
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Fig. 3. The dissipative voltage, for an electron gas sample, 
calculated from Eq. (3) as a function of E y . The parameters 
used refer to the experimental conditions of [2]. 

effective Euler equation, derived by Stone [7] , 

m* [v - [v x fl]] = e (E + v x B)- V |v| 2 + , 

(4) 

where v is the velocity field of the QHF, /1 is the local 
chemical potential containing all the interaction terms 
and is the fluid vorticity. Examining small pertuba- 
tions in the velocity and combining Eq. (4) with the 
continuity equation for the density of the QHF we find, 
in the limit \/fi — 0, 

hoj — Ae = — eE y lB(lBQ) + hojc- (5) 

For Ae = 0, this result is equivalent to Eq. (2) in the 
absence of interactions and corresponds exactly to the 
elastic inter-Landau level tunneling condition for the 
breakdown of the integer QHE [4] . Alternatively within 
this hydrodynamic model it corresponds to the condi- 
tion required to generate a vortex-antivortex pairs at 
zero energy. From Eq. (5), we can deduce that for a 
fluid described by Eq. (4) laminar flow becomes unsta- 
ble when ^ 

Ey = el B (l C B QY (6) 
To make a direct comparison with our earlier quan- 
tum mechanical calculation we need to evaluate the 
dissipative voltage drop along the sample due to the 
generation of these vortex-antivortex pairs from a sin- 
gle impurity. For a specific system we would have to 
rely on a numerical simulation of Eq. (4). However, 
we can make considerable progress by considering vor- 
tex shedding in classical fluid mechanics [8]. Consider 
an obstacle in the path of a fluid. At a low fluid ve- 
locity the flow around the obstacle is laminar. When 
the flow rate is increased, vortex-antivortex pairs are 
formed in the vicinity of the obstacle. However, a vor- 
tex street is not formed until the flow is fast enough 
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to free the vortex-antivortex pairs from the local flow 
field near the obstacle. In this steady state, each vortex- 
antivortex pair moves away from the obstacle at a ve- 
locity which is governed by the background fluid veloc- 
ity. This analogy suggests that the vortex-antivortex 
pair in a QHF moves away from the impurity at a veloc- 
ity given by E v /B. From classical hydrodynamics [8] 
it is also known that the distance between each vortex- 
antivortex (I) pair generated is approximately three 
times the separation between a single vortex and an- 
tivortex (d) (d/l = 0.28). Now consider two states for 
our fluid: firstly, the state where the fluid flow is lami- 
nar; then we know that the generation rate of vortex- 
antivortex pairs is zero. When Ae = it cost no en- 
ergy to generate excitations and the flow ceases to be- 
come laminar due to the formation of vortex-antivortex 
pairs. Hence from Eq. (6) we can evaluate E v /B. Divid- 
ing this velocity by the distance between each vortex- 
antivortex pair (d = 0.281b(IbQ)), to obtain the rate 
of generation of vortex-antivortex pairs, we obtain the 
voltage drop along the Hall bar as 

_ 0.28^ e 
Vx ~ e(l B Q) 2 ■ (V 

To estimate the dissipative voltage step height from 
this hydrodynamic model we have to evaluate Eq. (7). 
We can do this by referring back to the previous quan- 
tum mechanical calculation to give us an estimate of 
QIb- Taking the value for QIb = 1.9 for which the 
quantum calculations give the maximum value for V x 
we find, using Eq. (9) that AV X = 4.8mV for the NIST 
experiment [2]. 

In conclusion we have examined two models for the 
step-like breakdown of the integer QHE. The first is 
a full quantum calculation of the dissipative voltage 
drop along the Hall bar generated by scattering from a 
charged impurity. The second method uses an effective 
fluid model to describe the QHF from which we derive 
the elastic inter-Landau level tunneling condition [4]. 
From this condition and drawing analogies with clas- 
sical fluid dynamics we also obtain AV X . In both cases 
we compared our theory to the experiments of Lavine 
et al. [2] which were carried out on the U.S. resistance 
standard sample. In each case the agreement between 
theory and experiment is quite good. It should also be 
noted that we have also made comparisons of our mod- 
els with the step-like breakdown observed for a QHF 
of holes [9] ; again our theoretical results match exper- 



iments quite well [10] . 

The models which we have presented above are for 
a specific type of breakdown in which V x increases as a 
series of steps of regular height. However, many sam- 
ples do not exhibit such behaviour at breakdown. In 
such samples a single large increase in V x is observed at 
breakdown [5]. This behaviour is generally explained 
by the BSEH model [5] , which describes an avalanche 
process which excites electrons over an extensive area 
of the sample. The model presented above does not 
contradict this model; in fact it complements it. Our 
model which describes a local breakdown due to scat- 
tering from charged impurities, can be a precurser to 
this avalanche breakdown which spreads over much of 
the sample, i.e. in some samples it is possible for the 
conditions for local charged impurity-induced break- 
down to be met before the avalanche process takes over. 
The natural question then becomes what are these con- 
ditions and when can they be met. The condition for 
local charged impurity induced breakdown is sample 
specific; we require a locally high electric field in the 
region of a charged impurity to generate a stream of 
e-h pairs as compared to the total electric field across 
the sample reaching a critical value to trigger BSEH. 
We suggest that in the NIST samples some mechanism 
is inhibiting the avalanche, thus allowing us to observe 
the weaker local step-like breakdown. 
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